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ABSTRACT 

The three bright TeV blazars Mrk 421, Mrk 501 and PKS 2155-404 are highly 
variable in synchrotron X-ray emission. In particular, these sources may exhibit 
variable time lags between flux variations at different X-ray energy bands. How- 
ever, there are a number of issues that may significantly bias lag determinations. 
Edelson et al. (2001) recently proposed that the lags on timescales of hours, 
discovered by ASCA and BeppoSAX, could be an artifact of periodic gaps in the 
light curves introduced by the Earth occultation every ~ 1.6 hr. Using Monte 
Carlo simulations, in this paper we show that the lags over timescales of hours 
can not be the spurious result of periodic gaps, while periodic gaps indeed in- 
troduces uncertainty larger than what present in the evenly sampled data. The 
results also show that time lag estimates can be substantially improved by using 
evenly sampled light curves with large lag to bin-size ratio. Furthermore, we 
consider an XMM-Newton observation without interruptions and re-sample the 
light curves using the BeppoSAX observing windows, and then repeat the same 
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cross correlation function (CCF) analysis on both the real and fake data. The 
results also show that periodic gaps in the light curves do not significantly distort 
the CCF characters, and indeed the CCF peak ranges of the real and fake data 
overlap. Therefore, the lags discovered by ASCA and BeppoSAX are not due to 
periodic gaps in the light curves. 

Subject headings: BL Lacertae objects: general — galaxy: active — galaxy: 
nuclei — numerical methods — X-rays: galaxy 

1. Introduction 

One of the main advances from recent X-ray observations of blazars is the discovery 
of the energy-dependent time lags in the X-ray emission of the three bright TeV-emitting 
blazars Mrk 421, Mrk 501 and PKS 2155-304. The overall spectral energy distributions 
(SEDs) show that the synchrotron emission component from these sources peaks at high- 
energy (UV/soft X-ray) band. This indicates that the X-ray emission from these sources 
is the high energy tail of the synchrotron component, where the most violent variability is 
expected. The intcr-band time lag is one of the important variability parameters. The cross- 
correlation function (CCF) technique is the standard tool for lag determinations. ASCA and 
BeppoSAX discovered that in these sources the lower energy X-ray photons may lag or lead 
the higher energy ones, i.e., the so-called soft or hard lag, respectively. The signs and values 
of the lags may well depend on either energy or on the single flare analyzed. The typical 
lags (either soft or hard lag) range from ~ zero to 10^ s (see Zhang et al. 2002 for a review), 
and appear to be correlated to the flare duration: the shorter the flare duration, the smaller 
the lag (Zhang et al 2002; Brinkmann et al. 2003). Zhang (2002) also found evidence for a 
dependence of lags on timescalcs in Mrk 421: lags appear to be larger on longer timescales. 
The observed lags have been interpreted as evidence for the interplay of the acceleration and 
cooling timescales of relativistic electrons responsible for the observed X-rays taking place in 
the jets (e.g.. Kirk, Rieger & Mastichiadis 1998), and used to constrain the main parameters 
of the emitting region on the basis of variabihty (e.g., Zhang et al. 2002) rather than of 
SEDs, as commonly used. 

However, the latest lag searches with XMM-Newton have put the discovery of ASCA 
and BeppoSAX into question. Edelson et al. (2001) and Sembay et al. (2002) reported that 
there is no evidence for measurable inter-band lags in the X-rays, with upper limits of ~ 
0.3 (PKS 2155-304) and 0.08 hr (Mrk 421). It is commonly believed that uninterrupted 
time series, due to the highly eccentric ~ 48 hours orbit of XMM-Newton, allow to detect 
more reliably the lags. This has led these authors to propose that the lags measured with 
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ASCA and BeppoSAX could well be an artifact of the periodic interruptions by the Earth 
occultation related to the short ASCA and BeppoSAX oihital period (~ 1.6 hr). Maraschi 
et al. (2002) also reported no measurable lags with a full orbit XMM-Newton observation 
of PKS 2155-304. Brinkmann et al. (2003) re-analyzed all the available XMM-Newton 
observations of Mrk 421, working with light curves that have a small bin-size of 8 s, and 
found typical lags of ~ 5 min. 

In order to address the issue of the reliability and significance of the lags discovered by 
ASCA and BeppoSAX, in this paper wc perform Monte Carlo simulations to investigate the 
effects of periodically gapped time series on CCF lag determinations. A number of simula- 
tions using inverse Fourier transformation, which is model-dependent, have been carried out 
in the literature to evaluate the significance and uncertainty of the CCF lag determinations 
either for poorly sampled time series (e.g., Gaskell & Peterson 1987; Maoz & Netzer 1989; 
White & Peterson 1994) or for evenly sampled time series (e.g., Welsh 1999). The issue of 
poorly irregular sampling and noisy data was also previously investigated (see the summary 
in Kocn 1994). There are a number of issues that can significantly bias the reliability of CCF 
lag determinations. We refer to Welsh (1999) for a review on such issue, where the author 
also showed the effects of power spectral density (PSD) slope, duration, signal-to-noise ratio, 
de-trending and tapering of light curves. We note, however, that the previous simulations 
mainly dealt with the relationship between the UV/optical emission line flux and the contin- 
uum variations in Seyfert galaxies, which applied to the reverberation mapping (to measure 
the size of broad-line region) and typical ground-based irregular sampling patterns. In the 
blazar context, the simulations by Litchfield, Robson, & Hughes (1995) were performed on 
the radio sampling patterns. Our simulations are tailored to match the sampling charac- 
teristics of space-based observations introduced by low earth orbit satellites, which to our 
knowledge have not been explicitly addressed in astronomical literature. Another distinction 
is that wc discuss the time lags between the variations of the synchrotron emission in different 
X-ray bands. In particular, it is worth noting that there are obvious difficulties in accurately 
determining lags in the synchrotron X-ray emission of TeV blazars as (1) the lags are short 
compared to the bin-sizes of the available data; (2) the data are usually not equally sampled 
- in particular periodical interruptions of space-based observations are unavoidable for the 
low Earth orbit satellites; and (3) there are ambiguities in interpreting the complexities of 
the CCF results. 

In this work we also make use of the method introduced by Peterson et al. (1998), 
namely a model-independent Monte Carlo method to assess the uncertainties in the lag 

measurements obtained. This method does make use of real data, and is known as fiux 
redistribution/random subset selection (FR/RSS) that is based on the "bootstrap" method. 
Here we adopt it to study the effects of periodic gaps in lag determinations. To do so, we 
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use a real XMM-Newton observation without interruptions, and re-sample it with the typical 
BeppoSAX sampling windows. The same CCF analysis is then performed on both the real 
and fake time series. 

In §2 we conduct model-dependent simulations: §2.1 illustrates the assumptions for the 
simulations to be performed; and the results are presented in §2.2. The model- independent 
simulations are performed in §3. We discuss the significance of our results in §4. 

2. Monte Carlo Simulations 
2.1. Assumptions 

The inverse Fourier transformation from frequency to time domain is usually performed 
to simulate light curves by assuming a PSD model. We use the algorithm of Timmer & 
Konig (1995), which randomizes both the amplitude and the phase of a red noise process 
at each Fourier frequency. In exploring the full variety of possible light curves showing the 
same PSD, this algorithm is superior to the commonly used one that randomizes phases 
only (Benlloch et al. 2001; Uttley, M'^Hardy, & Papadakis 2002). In order to mimic a 
real situation as closely as possible, we normalize fake light curves on the basis of the real 
variabihty behavior of one observation of Mrk 421 obtained with BeppoSAX (Fossati et al. 
2000; Zhang 2002) as this observation has a good signal-to-noise ratio that will lay stress on 
the key point of this investigation. The assumptions and procedure of our simulations are 
as follows. 

A simple power-law PSD, i.e., P{f) oc /~" with slope a = 2.5 is assumed (Zhang 2002) 
to represent the red noise variability of the bright TeV blazars (Kataoka et al. 2001; Zhang 
et al. 1999; 2002), and from this the light curves are recovered on the basis of Monte Carlo 
technique. With one set of Gaussian distributed random numbers, we construct a fake pair 
of light curves that evenly sample 171 points with bin-size At = 512 s. This choice resembles 
the binned light curves of Mrk 421 used in Zhang (2002) to perform the CCF analysis. At 
the same time, in order to mimic the discovered time lags, we delay the phase of the second 
light curve of the fake pair on the basis of the relationship A0(/j) = 27r/jr(/i), where fi is 
the Fourier frequency (i = 1, A^/2; N — 171 is the total number of the evenly sampled 
light curve points), T(/j) is the time lag and A0(/j) the phase lag at fi. For simplicity, we 
assume that r is frequency-independent, i.e., r = constant (see, however, Zhang 2002 for the 
evidence of the dependence of r on frequency), and in turn A0(/j) = 27rr/j. An assumed 
(true) lag r therefore corresponds to the lag determined with the CCF methods. The aim of 
our simulations is to recover the assumed r by applying the CCF methods to the fake pair 
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with specific sampling windows. The two hght curves of the fake pair are then scaled to have 
the same mean and variance as the real 0.1-2 keV and 2-10 keV light curves of Mrk 421 
(Zhang 2002), respectively. In order to mimic photon counting (Poisson) white noise, the 
two fake light curves are further Gaussian randomly redistributed on the basis of the average 
errors on the real 0.1-2 keV and 2-10 keV light curves of Mrk 421. In order to simulate 
the periodic gaps of the light curves obtained with BeppoSAX, we re-sample the two fake 
hght curves by applying the real observing windows on the 0.1-2 keV and 2-10 keV hght 
curves of Mrk 421. Three CCF methods, i.e., interpolation cross correlation function (ICCF, 
White & Peterson 1994), Discrete Correlation Function (DCF, Edelson & Krolik 1998), and 
Fisher's z-transformed DCF (ZDCF, Alexander 1997) that is based on the DCF, are then 
used to cross-correlate the two fake light curves before and after applying the real observing 
windows. All CCFs are normalized by the mean and standard deviation of the two cross- 
correlated light curves using only the data points that actually contribute to the calculation 
of each lag (White & Peterson 1994), since the light curves in our cases arc not stationary. 
A simulation is deemed to have succeeded if rmax (the maximum value of the CCF) between 
the two fake light curves is significant at a level of confidence greater than 95% (Peterson 
et al. 1998). For each succeeded trial, we record the lags using three techniques to interpret 
the CCF results: (1) using the lag corresponding to rmax of the CCF, Tpeak; (2) computing 
the centroid of the CCF over time lags bracketing rmax, Tccnt ^ all the CCF points with r 
in excess of O.Srmax are used (Peterson et al. 1998); (3) fitting the CCF with a Gaussian 
function to find the location of the CCF peak, rgt. We repeat this procedure 2000 times to 
construct probability distributions of the lags (i.e., the cross-correlation peak distribution, 
CCPD; Maoz & Netzer 1989) for Tpeak, Tcent, and rgt. 



2.2. Results 

We simulate time lags covering the range discovered by ASCA and BeppoSAX in TeV 
blazars. Therefore, we arbitrarily assume r = 0, 300, 1400, 3000, 5400, 7100 s, respectively. 
We create fake light curves using the same sets of random number for each case except for 
different r. 

Figure 1 example, shows a pair of fake light curves that are evenly sampled, and 

the second light curve (solid circles) is delayed by r = 3000 s with respect to the first one 
(open circles). The lag is clearly visible by comparing both the peaks and the troughs of 
the two light curves. In this case, the ratio of r to At of the light curves is ~ 6. The ICCF, 
DCF and ZDCF of the two hght curves are shown in Figure Ib-d (solid hnes). We use Tpeak, 
Tcent, and Tfit to measure the lag, respectively. The results are reported in Table 1. One 
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can see that the true lag is properly recovered by the different CCF methods and different 
techniques used to interpret the CCFs. 

We then re-sample these two light curves using the real BeppoSAX observing windows 
as mentioned in §2.1. The resulted hght curves thus are affected by periodic gaps resembhng 
the light curves obtained with BeppoSAX. Note that the number of points (63) in the first re- 
sampled light curve is smaller than that in the second one (76) because on-board BeppoSAX 
the LEGS detector is less exposed than the MEGS one. The IGGF, DGF and ZDGF of the two 
light curves with periodic interruptions arc shown in Figure Ib-d (dotted line or open circles 
with error bars). The measured lags arc also reported in Tabic 1. In this case, both Tpcak and 
Tcent underestimate the true lag while rgt of the DGF and ZDGF still recover it, although 
with larger uncertainties. Figure Ib-d also show that the overall characteristics of the CCFs 
obtained from the evenly sampled light curves are almost identical to those obtained from 
the corresponding light curves with periodic gaps. However, due to the missing of a large 
number of data points in the latter case, the GGF peaks somewhat shift to smaller lags. 
Note that Tat is not calculated for the ICCF case because there are no estimates of the ICCF 
errors, and the DGF errors are overestimated with respect to the ZDGF ones. 

The statistical significance of the effects of periodic gaps on the GGF lag determina- 
tions is deduced from the probability distributions of a number of simulations. We show 
in Figure 2 (IGGF), Figure 3 (DGF), and Figure 4 (ZDGF) the GGPDs of the simulated 
lags for the case r = and 3000 s. It is worth noting that the GGPDs are almost always 
non-normal distributions, in particular for Tpcak and Tccnt- Therefore, we use the median and 
68% confidence level (with respect to the median) to statistically characterize the GGPDs. 
The statistical results for all the simulations performed are tabulated in Table 2. The first 
column is the assumed true lag, and columns 2-4 and 5-7 give the GGPD median and 68% 
confidence range for evenly sampled and periodically gapped light curves, respectively. As 
shown in Figures 2-4 and Table 2, the main results of our simulations can be summarized as 
follows: (1) the assumed true lags are recovered in all cases in terms of the GGPD medians 
with 68% confidence errors; (2) the main effect of periodic gaps is to broaden the GGPDs, 
thus to increase the uncertainty of the GGF lag determinations. Most importantly, periodic 
gaps do not produce artificial GGF lags. More specifically, from the simulations with true 
small lags (e.g, the cases of and 300 s lags), one can see that periodic interruptions in the 
light curves definitely do not produce spurious lags on timescales of hours; (3) in some cases, 
the three GGF methods and the three techniques used to interpret the GGF results do not 
give rise to completely consistent results. For example, Tpeak strongly depends on the lag 
steps used to calculate the GGF; the errors of the DGF are overestimated, thus producing 
the broadest GGPDs. 
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These results imply that (1) evenly sampled light curves with large lag to bin-size ratios 
give more reliable CCF lag determinations; (2) more importantly, the lags on timescales of 
hours discovered by ASCA and Beppo SAX can not be an artifact of periodic gaps in the light 
curves that have intrinsically small lag; (3) the only effect of periodic gaps is to introduce 
uncertainty and to increase the variance on the CCF lag determinations. 

3. A Specific Case: XMM-Newton Observations of PKS 2155-304 

The simulations presented in the previous section showed that periodical gaps in the 
light curves do not produce spurious lags between them, but only increase the uncertainty 
and variance in lag determinations. Uninterrupted data with high temporal resolution are 
available from XMM-Newton observations of Mrk 421 and PKS 2155-304. All of the CCF 
analysis, performed by Brinkmann et al. (2001; 2003), Edelson et al. (2001), Sembay et 
al. (2002), and Maraschi et al. (2002) showed that the inter-band lags between the soft 
and hard energy band are close to zero, with upper limits of about 1000 s. As we pointed 
out in the Introduction, these results led Edelson et al. (2001) to suggest that previous 
claims of time lags on time scales of hours might be an artifact of the periodic interruptions 
every ~1.6 hours due to the low-Earth orbits of satellites such as ASCA and BeppoSAX. 
Therefore, in addition to the simulations presented in the previous section arguing against 
the above suggestion, an important test to assess the role of periodic gaps is to consider 
the XMM-Newton data and re-sample them according to the ASCA or BeppoSAX ohserving 
windows, and then repeat the CCF analysis on the fake ASCA or BeppoSAX data. Whether 
or not inter-band lags would be detected from the fake data would be a strong argument in 
favor or against the claim by Edelson et al. (2001). 

In order to perform such a test, we take the first part of the XMM-Newton observation 
of PKS 2155-304 (Maraschi et al. 2002). The details of the data reduction and CCF analysis 
will be presented in Maraschi et al. (in preparation). For our purposes we just extracted 
the light curves in two energy bands, i.e., 0.7-1 and 1-2 keV. The hght curves (without 
interruptions) are shown in Figure 5a. First, we performed the CCF analysis on them: the 
results are shown in Figure 5b-d (solid line) for the ICCF, DCF and ZDCF, respectively. We 
then re-sampled the light curves with the typical BeppoSAX sampling windows, obtaining 
two light curves with periodic interruptions resembling real BeppoSAX observations. We 
performed the same CCF analysis on the fake light curves with periodic gaps just as we did 
on light curves without gaps. The results are also shown in Figure 5b-d (dotted line or open 
circles with error bars). In all cases the two light curves are highly correlated near zero lag, 
and the CCFs calculated from the real and fake data overlap near their peaks, suggesting 
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that the periodic interruptions do not change the CCF character near the peak. We also 
measured the lags with the three techniques. The results are tabulated in Tabic 3. Due 
to the complexities of the CCFs, the measured lags depend on the CCF methods and the 
techniques used to quantify the CCFs. However, the results show in general that the real lag 
may be close to zero, with 0.7-1 keV photons very marginally leading the 1-2 keV photons. 

Finally, we performed FR simulations (i.e., Gaussian randomly redistributing the light 
curves on the basis of the quoted errors). The resulting CCPDs are shown in Figure 6. 
We used the same statistical method as we did in § 2.2 to characterize the CCPDs. The 
statistical results are tabulated in Table 4. Within the 68% confidence errors, Tpcak gives lags 
consistent with zero, while Tcent and Tfit suggest positive lags of ~1000-2000 s, but at a low 
confidence (~ 2a). Note that the low confidence of the lag detections may be caused by the 
complexities of the variability behavior and the dependence of variations on the energy band 
considered. This also explains the fact that different methods and techniques may give rise 
to inconsistent results. In any case the comparison of the results obtained - using the same 
method and technique - from the light curves with periodic gaps and without gaps does not 
favor the suggestion that the periodic interruptions may produce spurious lags. 



4. Discussion and Conclusions 

We performed two sets of simulations to investigate the effects of space-based obser- 
vations with short orbital period (~ 1.6 hr) on the reliability and significance of CCF lag 
determinations, specifically for the energy-dependent variations of synchrotron X-ray emis- 
sion of TeV blazars. The first set of simulations (§2) make use of fake light curves generated 
with the Fourier transformation method. We investigated two main issues: (1) the effects 
of periodic data gaps in the light curves; and (2) the eff'ects of different lag to bin-size ra- 
tios of light curves. The simulations showed that evenly sampled fight curves indeed yield 
more reliable CCF lags, with smaller variance, than the periodically gapped light curves do. 
However, the CCPD analysis clearly showed that the light curves with periodic gaps still 
preserve the nature of the true lags (regardless of the values) even though they introduce 
larger lag variances than the light curves without gaps. Moreover, larger ratio of lag to light 
curve bin-size can improve the significance of CCF lag determination. The second set of sim- 
ulations is based on a real XMM-Newton observation without interruptions. We re-sampled 
it with the typical BeppoSAX sampling windows in order to study whether or not periodic 
interruptions may give rise to strong biases in CCF lag determinations. The complex nature 
of the variability in TeV blazars results in obvious difficulties when quantifying lags, and 
it is likely to be the source of discrepancies between the results quantified with different 
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techniques. However, the comparison of the results derived with the same CCF method 
and the same quantifying technique, shows that the hght curves with periodic interruptions 
do not produce spurious lags on timescales of hours if their intrinsic lag is indeed small. 
Therefore, our investigations argue against the proposal by Edelson et al. (2001) that the 
lags of about an hour discovered by ASCA and BeppoSAX are an artifact of periodic gaps 
introduced by low Earth orbit satellites. We thus conclude that the lags discovered by ASCA 
and BeppoSAX are most likely due to intrinsic variability properties of the sources, and not 
artificially produced from an intrinsic zero lag by periodic gaps. However, due to the com- 
plexities of variability that produces complicated CCF - irregularities and complexities of 
CCF could well be caused by different properties of variability on different timescales (Zhang 
2002) - it is not easy to quantify real lags with the CCF method and different CCF meth- 
ods and interpreting techniques may work for different cases. Our simulations also confirm 
that uninterrupted light curves with large lag to bin-size ratios can improve accuracy of lag 
determinations, in particular for small lags that require high sampling rates of light curves. 

Welsh (1999) showed that the reality of CCF lag determinations also depend (1) on the 
light curve auto-correlation function (ACF) sharpness (the sharper the ACFs, the narrower 
the CCF peak and the smaller the lag bias) and (2) on the ratio of the intrinsic lag to the 
duration of the light curves. The first dependence can be easily explained by recalling that a 
CCF is a convolution of two ACFs. The ACF sharpness is determined by the PSD steepness 
(the steeper the PSD, the broader the ACF). The X-ray PSDs of the three TeV blazars are 
steep (with slopes of ~ 2-3; Kataoka et al. 2001; Zhang et al. 1999; 2002), the ACF and 
the CCF peaks are therefore broad. De-trending hght curves might remove such bias, but 
it also removes the low- (Fourier) frequency variability of the sources. This can introduce 
serious errors into time series and needs to be done very carefully. The second dependence 
becomes important only for the lengths of light curves shorter than ~ 4 times the lag, which 
is never the case for the TcV blazars observed in the X-rays (typical ratios of the lengths of 
light curves to lags are 10-100). 

Finally, we stress that the real light curves of the bright TeV blazars are very complex 
(e.g., the relationship between light curves at different energies may not be represented by 
just one "fixed" lag). Such complexities definitely result in irregular CCF, e.g., the CCF 
peaks at zero lag but shows asymmetry, which makes the CCF methods less straightforward 
to lag determinations. The cross-spectral technique, a more complex tool used to determine 
the Fourier frequency-dependent lags, may have the advantage of avoiding such ambiguities 
at least over long timescales (see Zhang 2002 for details). However, because this method 
relies on the Fourier transformation, it is not applicable to unevenly spaced data. On the 
contrary, in time domain one can use the DCF to substitute the classical CCF when dealing 
with irregular data. 



-10- 



We thank the anonymous referee for constructive comments. This research has made 
use of the standard on-line archive provided by the BeppoSAX Science Data Center (SDC). 
This work is partly based on an observation with XMM-Newton, an ESA science mission. 
The Italian MIUR is acknowledged for financial support. The project is sponsored by the 
Scientific Research Foundation for the Returned Overseas Chinese Scholars, State Education 
Ministry. 

REFERENCES 

Alexander, T. 1997, in Astronomical Time Series, eds. D. Maoz, A. Sternberg and E.M. 
Leibowitz (Dordrecht:Kluwer), p. 163 

Benlloch, S., Wilms, J., Edelson, R., Yaqoob, T., & Staubert, R. 2001, ApJ, 562, L121 

Brinkmann, W., Sembay, S., Griffiths, R. G., Branduardi-Raymont, G., Gliozzi, M., Boiler, 
Th., Tiengo, A., Molendi, S., & Zane, S. 2001, A&A, 356, L162 

Brinkmann, W., Papadakis, I.E., den Herder, J.W.A., & Haberl, F. 2003, A&A, 402, 929 

Edelson, R. A., & Krolik, J. H. 1988, ApJ, 333, 646 

Edelson, R. A., Griffiths, G., Markowitz, A., Sembay, S., Turner, M. J. L., & Warwick, R. 
2001, ApJ, 554, 274 

Fossati, G., et al. 2000, ApJ, 541, 153 

Gaskell, CM., & Peterson, B.M. 1987, ApJS, 65, 1 

Kataoka, J., et al. 2001, ApJ, 560, 659 

Kirk, J., Rieger, F., & Mastichiadis, A. 1998, A&A, 333, 452 

Koen, C. 1994, MNRAS, 268, 690 

Litchfield, S.J., Robson, E.I., & Hughes, D.H. 1995, A&A, 300, 385 

Maraschi, L., et al., 2002, in New visions of the X-ray Universe in the XMM-Newton and 
Chandra Era, ed. F. Jansen, in press (astro-ph/0202418) 

Maoz, D., & Netzer, H. 1989, MNRAS, 236, 21 

Peterson, B. M., Wanders, I., Horne, K., CoUier, S., Alexander, T., Kaspi, S., & Maoz, D. 
1998, PASP, 110, 660 



Sembay, S., Edelson, R., Markowitz, A., Griffiths, R.G., & Turner, M.J.L. 2002, ApJ, 547, 
634 

Timmer, J., & Konig, M. 1995, A&A, 300, 707 

Uttley, P., M^Hardy I.M., & Papadakis, I.E. 2002, MNRAS, 332, 231 

Welsh, W.F. 1999, PASP, 111, 1347 

White, R.J., & Peterson, B.M. 1994, PASP, 106, 879 

Zhang, Y.H. 2002, MNRAS, 337, 609 

Zhang, Y.H., et al. 1999, ApJ, 527, 719 

Zhang, Y.H., et al. 2002, ApJ, 572, 762 



This preprint was prepared with the A AS IM^jX macros v5.2. 



- 12 - 



Table 1. CCF analysis results of one pair of fake light curves(s)^ 
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'^The true lag is assumed to be 3000 s. 



Table 2. CCPD analysis results of Monte Carlo simulation data (ks)^ 
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1Q+0.27 
-^-"-0.24 








7 

DCF 


-1 r|+0.40 

J^"-o.io 


7 


nq+0.48 

""^-0.38 






0.00 





00+°°° 

""-0.00 





nn+0.23 
""-0.23 





10t° 


18 

16 





00+°°° 

""-0.00 





nn+0.29 
""-0.07 





""-0.44 


0.30 





ri+0.00 

^-^-0.00 





9R+0.23 
^"-0.02 





30t° 


18 
17 





c-i+O.OG 
^■•■-0.51 





nr;+0.78 
"^-0.05 





90+0.47 
^^-0.45 


1.40 


1 


C4+0.00 
°*-0.00 


1 


qi +0.23 
'^^-0.02 


1 


381° 


28 
23 


1 


09+0.00 
"^-0.00 


1 


14+0.65 
-■-^-0.45 


1 


qfl+0.67 
^"-0.62 


3.00 


3 


07+0.00 
"' -0.00 


3 


"D-0.24 


2 


961° 


44 
38 


2 


09+0.00 
"^-0.00 


2 


47+0-74 
+ ' -0.41 


2 


Qn+0.85 

y^-0.79 


5.40 


5 


eq+0.00 
""^-0.51 


5 


qo+0.25 
''°-0.24 


5 


341° 


69 

64 


6 


1/1+0.00 
■■■^-1.54 


4 


qo+0.95 
'''-'-0.72 


5 


qr+0.78 
'-'^-0.62 


7.10 


7 


7+0.00 
-0.00 


7 


1 R+0.25 
^"-0.24 


7 


02l"o 


90 
81 


6 


(.(.+0.76 
""-0.00 


6 


01 +0-72 
^^-0.93 


7 


nQ+0.94 
"•J-0.82 
















ZDCF 












0.00 





""-0.00 





()Q+0.23 
""-0.23 





001° 


10 
10 





QQ + 0.00 

""-0.00 





Qf)+0.29 
""-0.07 





nQ+0.25 
"'-'-0.23 


0.30 





CI +0.00 
"^-^-0.00 





90+0.23 
^"-0.02 





311^ 


10 
10 





CI +0.00 
^■■■-0.51 





nc+0.78 
"^-0.05 





q/l+0.25 
'^^-0.25 


1.40 


1 


C/I+O.OO 
^^-0.00 


1 


01 +0.23 
"J ■'■-0.02 


1 


401° 


13 
12 


1 


no+o-oo 

"^-0.00 


1 


-1/1+0.65 
■"■^-0.45 


1 


40+0.37 
^"^-0.38 


3.00 


3 


07+0.00 
"' -0.00 


3 


Qp+0.24 
""-0.24 


3 


oolEJ 


18 

19 


2 


09+0.00 
"^-0.00 


2 


47+0-74 
+ ' -0.41 


2 


07+0.61 

' -0.54 


5.40 


5 


Ro+O.OO 
"•^-0.51 


5 


qo+0.25 
'->°-0.24 


5 


38115 


20 
32 


6 


1/1+0.00 
-■- + -1.54 


4 


qo+0.95 
^'->-0.72 


5 


40+0-48 
^-'-0.39 


7.10 


7 


1 7+0.00 
-0.00 


7 


1 f^+0.25 
■■■"-0.24 


7 


081S 


30 
41 


6 


f.f.+0.76 
""-0.00 


6 


qi+0-72 

''■^-0.93 


7 


-, q+0.49 

■'^"'-0.47 



"^The quoted values are the medians of the CCPDs, and the errors are 68% confi- 
dence range with respect to the medians. 
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Table 3. CCF analysis results of the XMM-Newton data (s) 







Real data 




Fake data 


"Tpeak 


Tcent 


Tfit 


'^peak 


Tcent 


Tfit 


ICCF 


300 


1904 




300 


-212 




DCF 





100 


695 ± 865 


1200 


707 


1538 ± 2547 


ZDCF 





100 


433 ± 179 


1200 


707 


879 ± 592 



Table 4. CCPD analysis results of the XMM-Newton data (s)"" 







Real data 






Fake data 






'7'peak 


Tcent 


Tfit 


Tpeak 


Tcent 


Tfit 


ICCF 

DCF 

ZDCF 


DUU_gQo 

r> + 0.600 

n+0.600 

^-600 


9021 +^^'' 
oni 1 +881 


1279^^^66 

719+418 
' ■'-^-358 


DUU_j200 

1200+6°° 
1 200+°°° 


1201+1°^"^ 

+ 1101 
' -1110 

1477+1101 

' -1110 


17671^7^7 

1 _g23 



'^The quoted values are the medians of the CCPDs, and the errors are 68% con- 
fidence range with respect to the medians. 
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-3x10* 3x10'* -SxlO'* 3x10* 

Time Lag (s) Time Lag (s) 

Fig. 1. — (a) An example of a fake light curve pair assuming r = 3000 s between them. For clarity, 
the count rates of the second light curve (solid circles) are lowered by 1.5. (b) the corresponding ICCF 
of (a). The solid line is obtained from evenly sampled light curves, and the dashed line is obtained from 
periodically gapped light curves rc-samplcd from (a) after applying the real BeppoSAX observing windows; 
(c) the corresponding DCF of (a). The solid line is obtained from evenly sampled light curves, and the open 
circles with error bars are obtained from the same periodically gapped light curves as in (b); (d) same as 
(c), but for the ZDCF. The dashed lines in (b), (c) and (d) indicate the true lag. 
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Fig. 2. — ICCF CCPDs for two different lags, (a) evenly sampled light curves with r = s; (b) periodically 

gapped light curves with r = s; (c) evenly sampled light curves with t = 3000 s; (d) periodically gapped 
light curves with r = 3000 s. The dotted line refers to Tpeak, and the solid line to Tcent- The vertical long 
dashed line indicates the true lag. 
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Fig. 3. — Same as Figure 2, but for DCF. The dotted line refers to Tpoak, the soUd line to Tcentj and the 
short dashed hne to rat- The vertical long dashed line indicates the true lag. 
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Fig. 4. — Same as Figure 2, but for ZDCF. The dotted line refers to Tpoak, the solid hne to Tcent, and the 
short dashed hne to rat- The vertical long dashed line indicates the true lag. 
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(a) Light Curves 
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Fig. 5. — (a) A specific case: light curves obtained with XMM-Newton for PKS 2155-304 (see text for 
details). The solid circles refer to the 0.7-1 keV energy band, and the open ones to the 1-2 keV band, (b) 

the corresponding ICCF of (a). The solid line is calculated from real data (i.e., evenly sampled light curves), 
and the dashed line from periodically gapped light curves re-sampled from (a) after applying the BeppoSAX 
sampling windows to (a); (c) the corresponding DCF of (a). The solid line is calculated from real data, and 
the open circles with error bars from the same periodically gapped light curves as used in (b); (d) same as 
(c), but for ZDCF. The dashed lines in (b), (c) and (d) indicate zero lag. 
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Fig. 6. — CCPDs obtained by Gaussian randomly redistributing the XMM-Newton light curves (Figure 5a) 
on the basis of the quoted errors. The upper panels are obtained from the real XMM-Newton light curves, 
and the lower panels from the periodically gapped light curves re-sampled from the real XMM-Newton light 
curves using the BeppoSAX sampling windows. The dotted line refers to Tpeak, the solid line to Tcent, and 
the short dashed line to TRf The vertical long dashed line indicates zero lag. 



